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-We  describ^an  algorithm  for  computing  the  skeleton  (medial-axis  sur¬ 
face)  of  an  object  defined  using  constructive  solid  geometry  (CSG)-  This 
surface  is  the  locus  of  all  points  in  the  object’s  interior  that  have  equal 
minimum  distance  from  at  least  two  distinct  parts  of  the  boundary.  The 
skeleton  can  be  used  in  blending,  motion  planning,  medical  tomography, 
computer  vision,  and  in  mesh  generation.  Wk  also  present  a  geometric 
analysis  of  Voronoi  surfaces  from  which  the  skeleton  is  composed,  f  ■  , , 

/  G 


1  Introduction 


A  common  design  paradigm  in  mechanical  engineering  is  to  create  complex 
mechanical  parts  from  primitive  shapes  by  regularized  Boolean  operations.  This 
paradigm  became  well-known  as  Constructive  Solid  Geometry  (CSG),  notably 

'Supported  in  part  by  NSF  Grants  CCR  86-19817  and  DMC  88-07550,  and  by  ON  It  Con¬ 
tract  N00014-86-K-0465. 


1 


•IWHmtmdN  arrJtTBfarr 

f-r  y.  jbtte 

i  *»  i'., ;  J •  v • :  ■  -r.  1/n  ilrnS  I «« 1 


(L 


with  the  work  by  Requicha  and  Voelcker,  e.gM  [13],  who  proposed  five  primitives: 
block,  sphere,  cylinder,  cone,  and  torus.  We  refer  to  all  objects  created  from 
the  named  primitives  as  CSG  objects.  They  provide  an  unambiguous  computer 
representation  of  manufacturable  objects  of  realistic  complexity  and  shape. 

We  consider  in  this  paper  the  formation  of  the  skeleton  of  CSG  objects, 
defined  precisely  below.  Skeletons  can  be  used  in  shape  representation  and  shape 
recognition  [1],  in  robot  motion  planning  [14],  and  for  the  purpose  of  defining 
certain  blending  surfaces  [3].  Moreover,  skeletons  can  be  used  to  decompose 
objects  as  a  first  step  in  mesh  generation  [10]. 

The  balance  of  this  section  defines  the  skeleton  and  its  constituent  elements, 
Voronoi  surfaces  and  Voronoi  curves.  It  also  gives  more  details  of  the  applica¬ 
tions  of  the  skeleton.  Section  2  develops  a  conceptual  approach  to  constructing 
the  skeleton  from  its  elements  and  explains  the  geometric  and  combinatorial 
properties  skeleton  construction  has  to  account  for.  A  general  procedure  feu 
constructing  Voronoi  surfaces  is  outlined  in  Section  3.  This  procedure  would  be 
used  in  situations  in  which  the  bounding  surfaces  of  the  CSG  object  are  compli¬ 
cated.  Simpler  surfaces  admit  special  constructions,  and  these  are  investigated 
in  Section  4,  from  a  geometric  perspective.  Section  5  concludes  the  paper. 


1.1  Skeleton  of  CSG  Objects 


The  locus  of  all  points  in  an  object  that  have  equal  minimum  distance  from  at 
least  two  bounding  faces  of  the  object  forms  the  interior  skeleton  of  I  lie  object. 
In  the  literature,  the  interior  skeleton  has  also  been  referred  to  as  the  medial-axis 
surface  [9]. 

Definition  1.1  (Blum)  The  interior  skeleton  of  a  3D  object  is  the  locus  of  the 
centers  of  all  its  maximal  inscribed  spheres.  The  interior  skeleton  of  a  2D  object 
is  the  locus  of  the  centers  of  all  its  maximal  inscribed  circles. 

The  dotted  lines  in  Figure  1  shows  the  skeletons  of  simple  2D  and  3D  objects. 
For  example,  the  skeleton  of  the  rectangular  block  shown  in  Fig  1(c)  consists  of 
13  planar  faces,  each  of  which  is  the  locus  of  ail  points  equidistant  to  a  pair  of 
bounding  faces  of  the  block.  Twelve  skeleton  faces  are  equidistant  to  a  pair  of 
adjacent  faces  of  the  block.  The  horizontal  skeleton  face  is  equidistant  to  the  lop 
and  the  bottom  faces  of  the  block.  Note  that  these  two  faces  are  not  adjacent. 
Therefore,  possibly  all  face  pairs  must  be  considered  when  constructing  the 
skeleton. 

The  faces  of  the  block’s  skeleton  can  be  thought  of  as  polygons  in  equidis¬ 
tant  planes.  To  construct  the  skeleton,  we  must  construct  the  individual  faces 
The  faces  lie  in  Voronoi  surfaces,  well-studied  in  computational  geometry.  We 
consider  them  next.  \ 
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Figure  1:  Skeletons  of  Sphere,  Rectangle,  Box 


1.2  Voronoi  Curves  and  Surfaces 

Voronoi  diagrams  are  widely  used  in  computational  geometry.  Given  a  set  S 
of  points  in  the  plane,  the  Voronoi  diagram  is  a  tesselation  of  t lie  plane  by 
polygons  such  that  every  point  p  of  S  is  contained  in  a  polygon  whose  inlerioi 
points  are  closer  to  p  than  to  any  other  point  of  5.  The  Voronoi  diagram  of 
a  finite  set  S  contains  closed  and  open  polygons,  where  an  open  polygon  is 
bounded  by  two  semi-infinite  edges.  Proximity  problems  such  as  closest  pair. 
Euclidean  minimum  spanning  tree,  triangulation,  etc.,  can  be  solved  eKicienlU 
when  using  Voronoi  diagrams  [12]. 

When  the  elements  of  S  are  simple  polygons  or  other  extended  two-dimensiona 
objects,  the  Voronoi  diagram  is  more  clearly  conceptualized  as  the  locus  of  cen¬ 
ters  of  all  maximal  circles  that  are  tangent  to  at  least  two  elements  of  S  from 
the  outside.  The  Voronoi  diagram  is  then  composed  of  curve  segments,  as  well 
as  straight  line  segments.  When  considering  extended  2D  objects,  a  similar 
process  wilhin  the  boundaries  of  each  object  yields  a  new  Voronoi  diagram;  one 
that  is  the  interior  skeleton  of  the  object.  We  therefore  speak  of  an  exterior  and 
an  interior  Voronoi  diagram.  The  exterior  Voronoi  diagram  can  be  thought  of 
as  the  interior  Voronoi  diagram  of  the  infinite  object  obtained  by  subtracting 
from  the  Euclidean  plane  the  interior  of  the  objects  in  S. 

If  the  elements  of  S  are  3D  objects,  we  obtain  a  higher-dimensional  Voronoi 
diagram  whose  faces  lie  on  surfaces.  As  in  the  two-dimensional  case,  wo  speak 
of  interior  and  exterior  Voronoi  diagrams.  Note  that  these  are  precisely  the 
interior  and  the  exterior  skeletons. 

When  the  set  5  contains  more  than  two  objects  or  one  object  whose  shape  is 
bounded  by  more  than  one  face,  the  skeleton  is  not  a  single  surface  but  consists 
of  fares.  Each  face  lies  on  a  certain  surface,  called  the  l  oronoi  surface.  For 
example,  the  Voronoi  surface  with  respect  to  a  plane  and  a  point  not  in  the 
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plane  is  evidently  a  paraboloid  of  revolution,  with  the  point  as  focus.  Voronoi 
surfaces  are  considered  in  Sections  3  and  4. 

1.3  Applications  of  Skeletons 

Object  skeletons  can  be  used  in  shape  representation  and  recognition  [  1  ].  Hr'iedy, 
the  skeleton  defines  a  decomposition  of  the  object  into  simpler  figures  that  can 
be  described  separately.  These  descriptions  are  then  combined  into  a  single, 
canonical  object  description  tiiat  is  independent  of  the  original  design  sequence. 
In  robot  motion  planning,  the  exterior  skeleton  of  a  set  S  of  3D  obstacles  can 
be  used  to  define  paths  that  maximize  the  distance  of  a  moving  robot  to  each 
obstacle;  [14]  uses  the  exterior  skeleton  to  find  collision  free  paths. 

In  computer-ciided  design  and  modelling,  variable-radius  ( rolling- bail )  blends 
can  be  defined  via  a  spine  curve  and  a  radius- variation  function.  The  exterior 
skeleton  surface  can  be  used  for  an  analytic  definition  of  the  spine  curve  [3]  when 
the  blend  is  a  fillet.  The  interior  skeleton  can  be  used  for  an  analytic  definition 
of  the  spine  in  the  case  of  rounds. 

In  [10]  it  is  argued  that  the  interior  skeleton  can  be  used  to  decompose 
objects  for  the  purpose  of  finite-element  mesh  generation  in  2D  domains  that 
have  been  defined  in  a  CSG  style  from  disks  and  rectangles.  Similarly,  we 
consider  decomposing  3D  CSG  objects  with  the  skeleton.  In  solid  mechanics, 
the  interior  skeleton  would  be  required.1  In  fluid  mechanics,  the  exterior  skeleton 
would  be  used,  e.g.,  when  studying  an  airfoil  in  turbulent  media. 

Despite  great  practical  demand  and  potential,  automatic  interfaces  of  sol  id 
modelers  with  analysis  codes  appear  to  be  severely  limited  [4].  The  ability  to 
generate  skeletons  automatically  from  a  solid  model  would  facilitate  construct¬ 
ing  more  sophisticated  interfaces  of  this  kind.  Algorithms  have  been  proposed 
for  the  two-dimensional  version  of  the  problem,  for  computing  interior  and  exte¬ 
rior  Voronoi  curves,  (11,  10].  Methods  for  constructing  simple  Voronoi  surfaces 
have  also  been  considered  (9,  3].  However,  complex  Voronoi  elements  (i.e.  com¬ 
plicated  curves  and  surfaces)  arise  when  skeletons  of  moderately  complex  CSG 
objects  are  considered,  and  there  seems  to  be  no  work  on  this  subject. 

In  this  paper,  we  take  a  deeper  look  at  such  skeleton  surfaces  by  investigating 
the  individual  Voronoi  elements  that  make  up  th*i  interior  skeleton  of  CSG 
objects.  It  is  shown  that  the  geometry  of  such  skeletons  are  inherently  complex. 
An  investigation  into  the  geometry  of  the  associated  Voronoi  surfaces  brings  to 
light  the  necessity  of  going  beyond  natural  quadrics  in  existing  solid  modellers. 

'For  an  automatic  method  to  generate  rectangular  grids  for  three-dimensional  polyhcdi  I 
domains  efficiently  see  [15]. 
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1.4  Preliminaries 


Our  usage  of  the  term  skeleton  refers  to  the  composite  of  Voronoi  curves  and 
surfaces  that  result  from  appropriate  pairs  of  bounding  faces  of  the  object. 
A  Voronoi  curve  is  defined  with  respect  to  a  pair  of  curves,  and  a  Voronoi 
surface  with  respect  to  a  pair  of  surfaces.  It  suffices  to  restrict  attention  to  the 
interior  skeleton,  since  the  exterior  skeleton  may  be  obtained  by  considering  the 
(regularized)  complement  of  the  object,  in  3-space. 

Mathematically,  a  Voronoi  surface  can  be  defined  as  follows.  Let  /  and  g  l>e 
two  surfaces.  Then  their  Voronoi  surface,  denoted  Vor (f,g)  is  the  locus  of  all 
points  equidistant  from  /  and  from  g: 

Vor(/,ff)  =  {p  6  R3  |  dj{p)  =  dg(p)} 

where  d/(p)  and  dg(p)  are  the  perpendicular  distances  of  point  p  from  f  and 
from  g .  respectively.  Similarly,  we  define  the  Voronoi  curve  of  curves  /  and  g 
as  the  locus  of  points  in  R2  that  are  equidistant  from  /  and  g.  If  /  or  g  are 
curves  in  3-space,  then  the  equidistant  points  form  a  surface,  again  called  the 
Voronoi  surface  of  /  and  g.  It  can  be  shown  that  Voronoi  curves  and  surfaces 
are  semi-algebraic  whenever  /  and  g  are. 

All  CSG  objects  are  bounded  by  faces  that  lie  on  planes,  natural  quadrics, 
and  tori.  When  the  entire  surfaces  are  considered,  e.g.,  the  entire  plane,  the 
infinite  cylinder,  the  infinite  double  cone,  etc.,  their  pairwise  Voronoi  surfaces 
will  usually  be  algebraic  surfaces.  However,  there  are  cases  in  which  a  Voronoi 
surface  is  only  semi-algebraic,  that  is,  it  is  a  subset  of  an  algebraic  surface.  This 
situation  arises  because  certain  Voronoi  surfaces  may  end  abruptly,  and  we  give 
an  example  of  this  phenomenon  later  in  Section  3. 


2  On  the  Construction  Skeletons 

The  skeleton  of  a  CSG  object  consists  of  Voronoi  elements.  These  elemenis 
are  obtained  by  considering  the  bounding  surfaces  of  the  object,  individually  or 
pairwise.  In  this  section,  we  illustrate  the  nature  of  the  skeleton  and  argue  that 
the  construction  of  the  skeleton  consists  of  trimming  Voronoi  surfaces  between 
pairs  of  boundary  elements  of  equal  or  mixed  dimensionality. 

2.1  Reflexive  Voronoi  Surfaces 

Reflexive  Voronoi  surfaces  are  lower-dimensional  Voronoi  surfaces  that  arise  in 
the  skeleton  of  certain  primitive  shapes.  Consider  a  sphere.  Its  skeleton  is 
just  a  point,  namely  the  center  of  the  sphere.  Since  such  a  sphere  could  be 
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Figure  2:  Block  with  blind  hole  and  the  partial  skeleton 

the  CSG  object  under  consideration,  it  is  necessary  to  identify  the  ivjlexivc 
Voronoi  surface  of  a  given  surface,  defined  to  consist  of  all  points  p  in  space 
with  the  property  that  there  exists  more  than  one  point  on  the  given  surface,  at 
a  minimum  distance  from  p.  It  is  evident  that  the  given  surface  must  necessarily 
be  non-planar. 

Examples  of  CSG  objects  whose  reflexive  Voronoi  surfaces  are  actually  the 
interior  skeletons  include  the  cylinder,  cone  and  torus.  For  the  cyclinder  and 
the  cone,  the  reflexive  Voronoi  surfaces  are  their  respective  axes.  In  the  case  of 
the  torus  it  is  the  spine  -  a  circle  that  is  the  locus  of  centers  of  the  generating 
spheres.  In  the  following,  we  will  refer  to  such  curves  as  Voronoi  surfaces  and 
not  Voronoi  curves,  because  they  arise  from  considering  surfaces.  Moreover,  for 
complicated  surfaces  they  need  not  degenerate  to  points  and  curves. 

2.2  A  Procedural  Construction  for  Skeletons 

Let  object  5  be  a  rectangular  block  with  a  blind  hole,  as  shown  in  Figure  2. 
Let  g  denote  the  cylinder  used  in  the  CSG  construction  and  /  the  plane  where 
the  hole  begins.  It  will  be  shown  later  that  the  skeleton  of  this  object  includes 
part  of  a  circular  cone  which  has  for  its  base  the  circle  of  intersection  of  g  and 
/.  See  also  Figure  2.  The  entire  cone  is  the  Voronoi  surface  of  the  plane  /  and 
the  (infinite)  cylinder  g.  In  principle,  the  construction  of  the  skeleton  of  this 
object  could  be  described  as  follows. 

1.  For  each  bounding  surface,  /,  of  5,  construct  the  reflexive  Voronoi  surface 
Vor(/,). 

2.  For  each  pair  of  bounding  surfaces,  /;,/*  of  5,  construct  the  Voronoi  sur¬ 
face  Vor  (/j,A). 
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Figure  3:  Odd  and  Even  Voronoi  Surface  of  two  Spheres 

3.  Trim  Vor(/,)  and  Vorf /_,,/*)  suitably  to  yield  the  interior  skeleton  of  S'. 

This  motivates  considering  the  Voronoi  surfaces  between  pairs  of  bounding  sur¬ 
faces  of  the  object. 

Note  that  this  approach  to  constructing  a  skeleton  is  too  simplistic.  Ii 
ignores  the  fact  that  the  presence  of  edges  and  vertices  may  cause  other  surlaco 
elements  to  appear  in  the  skeleton. 

2.3  Voronoi  Surface  Parity 

We  wish  to  determine  Voronoi  surfaces  for  pairs  of  adjacent  bounding  surfaces 
of  the  object.  Given  two  surfaces  /  and  g,  the  surface  Vor(/,  g)  consists  of  two 
components.  One  component  consists  of  points  that  lie  either  on  the  outside  of 
both  /  and  g  or  on  the  inside  of  them.  We  call  such  a  Voronoi  surface  the  even 
Voronoi  surface  of  f  and  g.  The  other  component  of  Vor(/,  g)  is  the  odd  Voionoi 
surface  of  /  and  g  and  consists  of  points  that  are  on  the  inside  of  /  and  outside 
of  g,  or  on  the  outside  of  /  and  the  inside  of  g.  Figure  3  shows  both  components 
for  two  spheres  of  equal  radius.  Here,  the  even  Voronoi  surface  is  a  plane  that 
contains  the  circle  in  which  the  two  spheres  intersect.  Points  on  this  plane  have 
equal  minimum  distance  to  either  sphere,  and  always  lie  on  the  outside  ol  both 
spheres  or  on  the  inside  of  both  of  them.  As  shown  later,  the  odd  offset  is  a 
prolate  spheroid,  i.e.,  an  ellipse  rotated  about  one  of  its  axes.  The  points  of  the 
spheroid  have  equal  minimum  distance  from  the  two  spheres,  but  arc  always 
outside  of  one  and  inside  the  other.  Which  component  of  the  Voronoi  surfaces 
is  needed  depends,  roughly  speaking,  on  whether  the  surfaces  bound  cavities  or 
protrusions.  More  precisely,  in  the  terminology  of  [2],  it  depends  on  whet  lieu 
the  primitives  contributing  the  two  faces  are  used  positively  or  negatively  iu  the 
CSG  tree. 


Figure  4:  Skeleton  at  the  end  of  a  blind  hole 

2.4  Mixed-Dimensional  Pairs 

Consider  the  bottom  of  the  blind  hole  of  Figure  2.  In  its  vicinity,  the  skeleton 
has  a  more  complicated  structure  as  illustrated  in  Figure  4.  Here,  the  skeleton  is 
composed  of  three  elements,  labeled  h i,  /i2,  and  /13  in  the  figure.  The  element  /i , 
is  part  of  a  cone,  the  Vor(/i,#)  of  cylinder  f\  bounding  the  hole  and  the  opposite 
face  plane  g  of  the  block.  The  surface  /13  is  part  of  a  plane,  the  Vor(iy,/2)  of  the 
plane  g  and  the  plane  /2  that  forms  the  bottom  of  the  hole.  The  intermediate 
surface,  /i2,  consists  of  points  that  have  equal  distance  from  the  plane  g  and  the 
edge  e  formed  by  the  intersection  of  f\  and  /2.  Thus,  h2  is  part  of  Vor(/i,c). 
We  will  see  later  that  Vor(/i,  is  a  parabolic  torus  —  a  surface  of  revolution 
obtained  by  rotating  a  parabc  iround  a  line  parallel  to  and  away  from  its  axis 
of  symmetry. 

This  shows  that  not  only  do  we  have  to  consider  Voronoi  surfaces  between 
pairs  of  bounding  surfaces,  but  also  between  surfaces  and  curves.  It  will  be  no 
surprise  that  we  also  need  to  determine  Voronoi  surfaces  for  curve  pairs,  for 
points  and  curves,  for  points  and  surfaces,  and  for  pairs  of  points.  Hence,  we 
have  to  augment  our  earlier  algorithm  with  a  determination  of  Voronoi  surface* 
of  such  pairs  of  boundary  elements,  of  lower  and  of  mixed  dimensionality. 

2.5  Trimming  Surfaces 

The  complexity  of  constructing  the  skeleton  of  CSG  objects  now  begins  to  show: 
Many  Voronoi  surfaces  have  to  be  constructed,  and  a  suitable  strategy  for  trim¬ 
ming  away  irrelevant  surface  areas  will  be  needed.  For  this  purpose,  we  should 
consider  where  the  “influence  region”  of  a  face  ends  on  an  associated  Voronoi 
surface.  This  is  especially  important  in  view  of  the  fact  that  since  adjacent 
components  of  a  skeleton  may  meet  with  tangent  continuity,  a  determination 
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Figure  5:  Trimming  Surface  at  Edge  e 

of  the  separating  edges  as  intersection  curves  of  the  associated  Voronoi  surfaces 
may  be  pose  delicate  numerical  problems. 

Consider  a  point  p  in  space.  We  ask  for  a  point  q  on  a  given  face  f  of  the 
CSG  object  that  lies  closest  to  p.  If  q  is  interior  to  /,  then  the  line  pq  must  be 
normal  to  /  at  q.  This  motivates  trimming  the  Voronoi  surface  associated  with 
f ,  by  a  surface  of  normals  at  the  boundary  of  /.  This  surface  is  a  ruled  surface 
generated  by  all  the  normals  to  /,  along  its  boundary. 

In  simple  situations  we  expect  that  the  trimming  curves  do  not  exhibit  com¬ 
plicated  singularities,  and  that  they  delimit  a  simply  connected  area  that  con¬ 
stitutes  a  face  of  the  skeleton.  Note,  however,  that  the  trimming  surface  could 
have  a  nontrivial  geometry.  For  example,  for  a  face  on  the  cylinder  /  bounded 
by  the  elliptic  edge  e,  the  trimming  surface  is  the  ruled  surface  /e  shown  in 
Figure  5. 

Similar  trimming  methods  can  be  formulated  for  edges:  Trim  the  associated 
Voronoi  surface  by  the  normal  planes  at  the  vertices  of  the  edges,  and  by  the 
trimming  surfaces  of  the  adjacent  faces.  In  summary,  we  construct  the  Voronoi 
surfaces  on  which  the  faces  of  the  skeleton  are  situated,  by  considering  elements 
of  the  boundary  individually,  and  in  pairs.  Then  Voronoi  surfaces  are  trimmed 
by  the  trimming  surfaces  which  are  generated  by  normals  to  the  surfaces  and  to 
the  edges.  The  skeleton  is  then  assembled  from  the  patches  of  Voronoi  surfaces 
that  remain. 

2.6  Algebraic  Complexity 

Aside  from  the  combinatorial  complexity  of  properly  trimming  and  combining 
patches  of  Voronoi  surfaces,  t lie  Voronoi  surfaces  themselves  can  be  algebraically 
quite  complex.  In  response,  we  present  in  the  next  section  a  genend  method 
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that  formulates  the  surfaces  in  a  fairly  simple  way,  as  a  set  of  equations  using 
auxiliary  variables.  These  equations  can  be  formulated  easily,  but  their  analysis 
and  interrogation  is  not  so  simple  and  is  the  subject  of  ongoing  research:  (0. 
7].  So,  we  give  a  simplified  method  for  deriving  implicit  equations  of  Voronot 
surfaces  of  CSG  surface  pairs.  In  addition,  in  Section  4,  we  discuss  the  geometry 
of  Voronoi  surfaces  for  CSG.  The  geometric  properties  derived  in  that  section 
provide  shortcuts  to  the  construction  and  analysis  of  the  elements  of  the  interior 
skeleton. 


3  Methods  for  Deriving  Voronoi  Surfaces 


Recall  that  the  Voronoi  surface,  Vor(/,  g),  is  the  locus  of  all  points  equidistant 
from  surfaces  /  and  g.  Given  /  and  g ,  implicitly  or  parametrically,  there  is  a 
general  method  for  formulating  a  mathematical  description  of  Voronoi  curves 
and  surfaces,  as  a  set  of  equations  using  auxiliary  variables.  Although  we  use 
this  general  method  for  algebraic  curves  and  surfaces  in  this  paper,  we  note  that 
it  can  be  used  for  all  curves  and  surfaces  that  are  differentiable.  We  recall  the 
method  from  [3,  5,  6,  7]. 

3.1  A  General  Procedure 

Let  /  and  g  denote  a  pair  of  algebraic  surfaces.  Consider  two  spheres  Sj  and 
Sg,  each  of  radius  d,  with  centers  on  /  and  g,  respectively.  With  d  the  offset 
distance,  the  ofTset  surfaces  of  /  and  g  are  given  by  the  envelopes  generated  by 
Sj  and  Sg  as  they  move  on  /  and  g.  This  is  equivalent  to  displacing,  in  the 
normal  direction,  the  centres  (uj,U2,U3)  of  5/  and  (vi,i>2,U3)  of  Sg ,  as  they 
move  on  /  and  g,  by  the  radius  of  the  spheres.  Denoting  the  pair  of  linearly 
independent  tangent  directions  to  /  at  any  point  by  tj  and  t2  and  .  imilai  ly  to 
g  by  t\  and  t2  we  formulate  the  following  eight  equations  in  ten  variables 

/:  f(uuu2,u3)  =  0 

9 ■  ff(vi,V2,u3)  =  0 

5/:  (x  -  ui)2  +  (y  -  u2)2  +  [z  ~  u3)2  -  d2  =  0 

Sg  :  (x  -  iq)2  +  (y  -  t>2)2  +  (*  -  v3)2  -  d.2  =  0 

(VS/  t,)  =  0 
(VS/-t2)  =  0 
(VSfft',)  =  0 
(VSj-t'j)  =  0 
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By  eliminating  the  variables  {u!,u2,  ti3,t>|,v2,u3,d}  the  implicit  equation  of  the 
Voronoi  surface  Vor (f,g)  associated  with  /  and  g  can  be  obtained,  at  least  in 
principle. 

An  example  helps  to  illustrate  the  above  procedure.  Let  /  be  a  cylinder 
of  unit  radius,  parallel  to  the  2-axis  with  the  line  (z  =  0 ,y  =  2)  as  its  axis. 
Let  g  be  another  cylinder  of  unit  radius,  parallel  to  the  z-axis  with  the  line 
(z  =  0,y  =  -2)  as  its  axis: 

/:  z2  +  (y-  2)1  -  1  =  0 

g  :  z2  +  (y  +  2)2  -  1  =  0 

Let  d  be  the  radius  of  the  spheres  Sj  and  Sg.  At  points  (itj,  112,143)  and 
( fi ,  ,  U3)  on  /  and  g  respectively,  we  have  the  equations 

f  ■  U?  +  (ti2  -  2)2  -1  =  0 

5/:  (z  -  «i  )2  +  (y  -  u2)2  +  (z  -  ii3)2  -  d2  =  I) 

9  •  u3  +  (v2  +  2)2  -  1  =  0 

Sg  :  (z  —  tq)2  +  (y  —  t>2)2  +  (*  -  u3)2  ~  d2  =  0 

The  gradients  V /  and  Vg  are  given  by  [2ui,2u2  -  4,0]  and  [0,2u2  +  4,2v3]. 
The  two  independent  tangent  vectors  at  (ui,u2,u3)  on  /,  are 

t,  =  (0,0,1)  t2  =  (2  -  u2,uj,0) 

At  (iq,u2,u3)  on  g  they  are 

t',  =  (1,0,0)  tj  =  (0,  —  v3,  u2  +  2) 

Thus,  the  final  system  of  eight  equations  is  given  by 

+  (u2  -  2)2  -  1  =  0 

(z  -  ui)2  +  (y  -  u2)2  -f  (z  -  u3)2  -  d2  =  0 

vl  +  (v2  +  2)2  -  1  =  0 

(x  -  Vi)2  +  (y  -  v-i)2  +  (z  -  V3)7  -  d2  =  0 

z  —  u3  =  0 

(z  -  Ui)(2  -  u2)  +  (y  -  U2)U!  =  0 

x  -  v\  =  0 

—{y  -  V?)V3  +  (z  -  v3)(t>2  +  2)  =  0 
After  elimination,  we  obtain  the  equation  of  the  Voronoi  surface. 
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Surface 

Equation 

r-Offset 

Equation 

Plane 

X 

Plane 

x  -  d 

Sphere 

x2  4-  y2  +  z2  -  r2 

Sphere 

x2  4-  m2  +  z2  —  ( r  4-  d)2 

Cylinder 

i 2  y2  —  r2 

Cylinder 

x2  -f  y 1  -  (  r  +  il)2 

Cone 

x2  +  y2  -  z 2  tan2(7) 

Cone 

x2  +  y2  -  {z  +  .//  si n( 7 ) )2  taii2( -  ) 

Torus 

(x2  +  y2  +  z2  -  r2)2 

Torus 

(x2  +  y2  +  z2  -  (  r  +  d)2)2 

+2c2(x2  -  (x2  +  y2  4-  r2))  +  c4 

+  2 c2(?2  -  (x2  +  y2  +  (r  +  d)2))  +  r1 

Table  1:  d- Offsets  of  CSG  Surfaces 

Elimination  of  the  additional  variables  in  an  attempt  to  derive  an  implicit 
equation  for  the  Voronoi  surface  is  a  task  well  beyond  the  capacity  of  cuno'it 
elimination  procedures  known  to  us;  see,  e.g.,  [7j.  Instead,  we  use  the  equations 
formulated  by  the  general  method  only  for  numerical  computations  that  analyze 
the  surface  locally,  as  in,  e.g.,  (5). 

We  mention  that  this  method  may  generate  extraneous  solutions.  Some 
of  these  may  be  eliminated  by  adding  additional  equations;  see  [6],  Whether 
such  additions  succeed  in  all  cases  is  not  known  at  this  time.  Note  also  that  the 
general  method,  as  described,  cannot  distinguish  between  odd  and  even  Voronoi 
surfaces.  In  consequence,  the  implicit  equation  derived  from  the  system  ( 1 )  must 
be  the  product  of  the  two  surfaces.  This  fact  would  be  unpleasant  for  a  global 
surface  analysis.  For  a  local  surface  analysis,  however,  the  additional  surface 
components  present  no  serious  obstacle. 

3.2  Simplifications  for  Special  Surfaces 

The  explicit  derivation  of  the  Voronoi  surface,  by  the  procedure  just  described, 
requires  excessive  symbolic  computations.  Considerable  simplifications  are  pos¬ 
sible  for  specific  classes  of  surfaces,  including  planes,  natural  quadrics,  and  tori. 
For  example,  if  one  of  the  primary  surfaces  is  a  plane  and  the  other  a  cylinder 
or  a  cone,  it  is  easy  to  see  that  the  associated  Voronoi  surface  is  a  cone  [3j.  c;ucli 
simplifications  are  based  on  the  idea  to  replace  the  four  equations  describing  the 
d-ofTset  of  a  surface  with  a  single  equation,  where  possible,  thereby  eliminating 
up  to  six  variables  outright.  Note  that  d-offsets  are  explicitly  known  for  planes, 
natural  quadrics,  and  the  torus.  They  are  summarized  in  Table  l.  Figures  <> 
through  8  explain  the  geometric  meaning  of  the  constants  in  the  equations.  Flic 
ofTset  equations  are  readily  verified. 

As  an  example,  we  derive  the  even  Voronoi  surface  of  th"  two  non- intersecting 
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Voronoj  Surface  of  two  Cylinders 


•v'  v 

*r 


S 


(Hyperbolic  parapoloicj) 


cylinders  discussed  before.  The  equations  for  the  respective  d-oflsets  are  then: 

F:  x2  +  (y  —  2)2  —  (1  -f  d)2  =  0 

G:  z2 +  (y  +  2)2 -(l  +  d)2  =  0 

Here,  d  is  eliminated  in  h  =  F  -  G,  which  is  the  sought  equation  of  the  even 
Voronoi  surface  of  the  two  cylinders: 

h  :  x2  -  8y  —  z2  =0 

It  is  a  hyperbolic  paraboloid,  as  shown  in  Plate-I. 

In  general,  d  has  to  be  eliminated  in  two  steps.  As  an  example,  consider  the 
odd  Voronoi  surface  of  the  two  cylinders.  It  is  described  by 

F  :  x2 +  (y-2)2 -(l  +  d)2  =  0 

G  :  z2  +  (y  +  2)2  —  (1  —  d)2  =  0 

We  obtain 

F  -  G  =  x2  —  8y  —  z2  —  4d 

from  which  we  determine  d  as  a  function  of  x,  y  and  z : 

d  =  (x2  -  8y  -  z2)/ 4 

Substitution  into  F  then  yields  the  equation  of  the  odd  Voronoi  surface  as 

h  =  x2  +  (y  -  2)2  -  (l  +  (x2  -  8y  -  z2)/ 4)2 

This  is  a  quartic  surface.  Since  the  cylinders  do  not  intersect,  physically  there 
cannot  exist  a  Voronoi  surface  with  an  odd  component.  Algebraically,  wo  might 
say  that  the  odd  Voronoi  surface  does  not  contain  any  real  points,  if  the  cylin¬ 
ders  have  a  nonempty  intersection,  the  odd  Voronoi  surface  will  contain  real 
points.  Such  an  example  is  shown  in  Plate-II. 

We  mentioned  in  Section  1  that  the  Voronoi  surface  of  two  algebraic  surfaces 
may  be  semi-algebraic.  This  is  the  case  with  the  odd  Voronoi  surface  shown  in 
Plate-II.  Here,  the  algebraic  surface  has  four  conical  singularities,  two  of  which 
are  visible.  The  central  part  of  the  surface,  in  shape  roughly  tike  a  pillow,  is 
the  odd  Voronoi  surface  and  ends  in  the  four  conical  points  that  lie  on  the 
two  cylinder  axes.  The  containing  algebraic  surface  continues  beyond  these 
singularities.2  Its  other  points  do  not  belong  to  the  Voronoi  surface,  since  the 
points  of  equal  distance  on  the  two  cylinders  minimize  the  distance  locally  only, 
not  globally. 

Hn  the  picture  the  four  conical  surface  extensions  have  been  clipped  almost  entiiclv. 
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4  Voronoi  Surfaces  between  Surface  Pairs 


We  now  consider  the  geometry  of  Voronoi  surfaces  that  result  between  pairs 
of  bounding  surfaces.  We  restrict  the  bounding  surfaces  to  be  planes,  natural 
quadrics,  and  tori.  Many  of  the  Voronoi  surfaces  that  arise  can  be  analyzed  by 
purely  geometric  argument,  and  we  do  this  where  appropriate.  Note  tiiat  we 
consider  the  bounding  surfaces  in  their  entirety.  The  Voronoi  surface  of  such 
a  pair,  as  pointed  out  earlier,  is  not  always  a  complete  algebraic  surface,  but 
is  always  part  of  an  algebraic  surface.  We  do  not  address  how  to  trim  Voronoi 
surfaces. 

First,  we  recall  the  relevant  geometric  facts  about  Voronoi  curves  lor  points, 
circles,  and  lines.  These  results  are  used  later  to  analyze  the  Voronoi  surlncos. 
In  the  following,  we  refer  to  planes,  natural  quadrics,  and  tori,  collectively  as 
CSC  surfaces. 

4.1  Voronoi  Curves  of  Points,  Circles,  and  Lines 
Proposition  4.1 

The  Voronoi  curve  of  a  straight  line  and  a  distinct  point  is  a  parabola.  If  the 
point  lies  on  the  line,  the  Voronoi  curve  is  the  perpendicular  to  the  line  at  that 
point. 

Proposition  4.2 

The  odd  and  even  Voronoi  curves  of  a  circle  C  and  a  straight  line  L  are  parabolas. 
Each  has  the  center  of  C  as  its  focus.  If  the  circle  is  tangent  to  L,  then  one  of 
the  parabolas  degenerates  into  a  half  line. 

Proof 

Let  the  radius  of  the  circle  C  be  r,  and  consider  first  the  case  where  L  intersects 
C.  Without  loss  of  generality,  we  orient  L  such  that  the  center  of  C  lies  to  (he 
outside  of  I.  or  on  L.  Consider  the  component  hi  of  the  even  Voronoi  curve  that 
lies  outside  both  L  and  C .  If  p  is  a  point  on  h\  at  a  distance  d\  from  the  circle 
C  and  the  line  L ,  then  p  is  at  distance  di  4-  r  from  the  center  of  the  circle.  Hence 
it  is  also  at  distance  dj  +  r  from  a  line  L\,  parallel  to  and  at  a  distance  r  I torn  L. 
on  the  inside  of  L.  See  also  Figure  9  left.  By  Proposition  4.1,  the  component  h\ 
of  the  even  Voronoi  curve  lies  on  a  parabola  P\.  Now  consider  the  component 
/i2  of  the  even  Voronoi  curve  that  lies  inside  C  and  inside  L.  If  p  is  a  point  on 
h2,  at  a  distance  d2  from  C  and  L ,  it  is  also  at  a  distance  r  -  d2  from  the  center 
of  C.  Hence  it  is  also  at  a  distance  r  -  d2  from  the  line  L i,  and  hence  is  also  on 
the  parabola  P\.  Note  that  the  parabola  passes  through  the  intersection  points 
of  L  and  C,  has  the  center  of  C  as  its  focus  and  the  perpendicular  to  L  passing 
through  the  center  of  C  as  its  axis. 
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Figure  9:  Voronoi  Curves  of  Line  and  Circle 

By  a  similar  argument,  the  odd  Voronoi  curve  can  be  shown  to  be  another 
parabola,  also  passing  through  the  intersection  points  of  L  and  C .  with  the  same 
axis,  but  an  opposite  orientation. 

Now  assume  that  L  is  tangent  to  C.  Then  the  component  /i2  of  the  even 
Voronoi  curve  no  longer  exists,  and  h\  and  the  point  of  tangencv  cover  the  entire 
parabola.  Moreover,  the  odd  Voronoi  curve  is  the  half  line  on  the  perpendicular 
to  L  at  the  point  of  tangency,  that  lies  on  the  inside  of  L.  See  also  Figure  9 
right.  When  L  and  C  do  not  intersect,  the  odd  Voronoi  curve  does  not  exist.  □ 

Proposition  4.3 

The  Voronoi  curve  of  two  nonintersecting  circles  is  a  hyperbola  if  the  circles 
have  different  radii,  and  a  straight  line  if  the  circles  are  of  equal  radius. 

Proof 

Let  C\  and  C2  be  the  two  circles,  and  consider  a  point  p  exterior  to  both,  on 
the  Voroni  curve.  If  d  is  the  distance  from  either  circle,  then  the  distance  to 
the  two  centers  is  d  4-  ri  and  d  +  r2,  respectively,  where  r,  and  r2  are  the  radii. 
Hence  the  difference  of  the  two  distances  is  constant.  If  the  radii  differ,  then  t  lie 
constant  is  not  zero,  and  hence  the  Voronoi  curve  is  a  hyperbola.  Otherwise, 
the  difference  is  zero,  i.e.,  the  Voronoi  curve  is  a  straight  line.  □ 

Proposition  4.4 

The  odd  Voronoi  curve  of  two  intersecting  circles  is  an  ellipse  and  the  even 
Voronoi  curve  a  hyperbola.  If  the  circles  have  equal  radii,  the  hyperbola  degen¬ 
erates  into  a  straight  line. 

Proof 

Let  Ci  and  C2  be  the  two  circles,  and  consider  a  point  p  exterior  to  both,  on  the 
even  Voronoi  curve.  If  d  is  the  distance  from  each  circle,  then  the  distance  to 
the  two  centers  is  d  -f  ri  and  d  4-  r2,  respectively,  where  ri  and  r2  are  the  radii. 
Hence  the  difference  of  the  two  distances  is  constant,  so  the  even  Voronoi  curve 


1G 


/ 

/  e*n. 

/ 


Figure  iO:  Voronoi  Curves  of  two  Circles 

is  a  hyperbola  if  r\  ^  r-i  and  a  straight  line  if  r\  =  r2.  Now  consider  a  point  (/ 
on  the  odd  Voronoi  curve,  with  distance  d  from  both  circles,  as  shown  in  Figure 
10.  The  distances  from  the  two  centers  are  therefore  rj  —  d  and  r2  +  </.  so  the 
sum  of  the  distances  is  constant.  Hence  the  odd  Voronoi  curve  is  an  ellipse  □ 

4.2  Voronoi  Surfaces  of  CSG  Surface  Pairs 

We  classify  the  types  of  Voronoi  surfaces  for  CSG  surface  pairs.  Although  the 
algebraic  degree  of  a  Voronoi  surface  can  be  as  high  as  16,  for  two  tori  in  general 
position,  there  are  many  cases  in  which  the  Voronoi  surface  is  much  simpler. 
Moreover,  there  are  many  geometric  properties  of  interest  to  be  reported. 

We  think  of  a  Voronoi  surface  as  the  projection  into  (x,y,z)- space  of  the 
intersection  of  two  hypersurfaces  in  (x,y,  z,r)- space.  Each  hypersurface  is  a 
generic  r-ofTset.  The  geometry  of  the  Voronoi  surface  may  depend  on  whether 
two  positive  r-o(Tsets  are  intersected,  or  whether  a  positive  and  a  negative 
r-ofTset  are  intersected.  We  speak  of  an  even  Voronoi  surface  in  the  posi- 
tive/positive  and  in  the  negative/negative  case.  In  the  positive/negative  and 
the  negative/positive  case,  we  speak  of  an  odd  Voronoi  surface. 

We  classify  first  those  Voronoi  surfaces  in  which  one  member  of  the  pair  is 
a  plane.  The  following  is  elementary: 

Proposition  4.5 

Consider  the  Voronoi  surface  Vor (f,g),  where  /  is  a  plane  and  (j  is  a  plane, 
sphere,  cylinder,  cone,  or  a  torus.  Then  the  degree  of  Vor (/,</)  docs  not  exceed 
the  degree  of  g. 

Proof 

Since  the  r-o(Tset  of  /  is  a  plane,  we  may  write  r  =  h(x,y,z),  where  h  is  linear. 
Consider  the  r-ofTset  of  g ,  where  g  is  a  CSG  surface.  Any  term  xy}  zkrm  in  the 


17 


✓ 


r-otfsct  equation  satisfies  i  +  j  +  k  +  m<  2  for  the  sphere,  cylinder,  and  cone, 
and  satisfies  i-\-j  +  k  +  m<4  for  the  torus.  Substitution  into  the  equation  for 
g  therefore  results  in  an  equation  of  degree  at  most  the  degree  of  g.  □ 

Note  that,  when  /  is  a  plane,  the  geometry  of  the  odd  and  the  even  Voro.ioi 
surface  do  not  differ.  This  is  analogous  to  Proposition  4.2.  The  individual 
Voronoi  surfaces  are  characterized  now: 

Theorem  4.6 

Let  h  —  Vor(/, g),  where  /  is  a  plane. 

(i)  If  g  is  a  plane,  then  so  is  h. 

(ii)  If  g  is  a  sphere,  then  h  is  a  paraboloid  of  revolution. 

(iii)  Let  g  be  a  cylinder.  If  the  axis  of  g  is  parallel  to  or  lies  in  /,  then  h  is  a 
parabolic  cylinder.  If  the  axis  of  g  intersects  /  in  a  point,  then  Ii  is  a  cone 
whose  baseline  is  the  elliptic  or  circular  intersection  of  g  and  /. 

(iv)  If  g  is  a  cone,  then  so  is  /.  If  /  contains  the  axis  of  g,  then  Ii  has  a 
hyperbola  as  base.  If  /  does  not  contain  the  axis  of  g ,  nor  its  vertex,  t  lien 
h  has  the  intersection  of  /  and  g  as  its  base,  as  base. 

(v)  Let  g  be  a  torus.  If  /  is  perpendicular  to  the  axis  of  the  torus,  then  h  is 
a  surface  obtained  by  rotating  a  parabola  about  the  axis. 

Proof 

Case  (i)  is  trivial.  For  Case  (ii)  observe  that  Vor (f,g)  must  be  a  surface  of 
revolution  about  the  normal  to  /  passing  through  the  center  of  g.  It  follows 
from  Proposition  4.2  that  h  is  a  paraboloid  of  revolution. 

Case  (iii):  Let  the  axis  of  cylinder  g  be  parallel  to  the  plane  /,  and  consider  a 
plane  P  perpendicular  to  the  axis  of  g.  P  intersects  /  in  a  line  L  and  y  in  a 
circle  C.  Because  P  is  also  perpendicular  to  the  plane  /,  the  points  on  }  and 
g  that  are  nearest  to  a  given  point  p  in  the  plane  P  must  also  lie  on  the  plane 
P.  So,  P  intersects  the  Voronoi  surface  of  /  and  g  in  the  Voronoi  curve  of  the 
circle  C  and  the  line  L ,  by  Proposition  4.2  a  parabola.  Since  the  position  and 
size  of  L  and  C  are  independent  of  the  position  of  P,  therefore,  Vor (/,</)  is  a 
parabolic  cylinder. 

Now  assume  that  the  axis  of  g  intersects  /.  Let  C  be  the  curve  of  intersection 
between  /  and  g.  Clearly,  C  is  either  a  circle  or  an  ellipse  and  lies  on  the  Voronoi 
surface  h.  Consider  the  plane  P  through  the  axis  of  g  and  perpendicular  to  /. 
P  intersects  /  in  a  line  L\  and  g  in  the  parallel  lines  Li  and  L -j.  Since  the 
closest  approach  on  /  and  g ,  to  any  given  point  p  on  P  must  lie  on  those  lines, 
P  intersects  h  in  two  lines  H\  and  Hi  that  bisect  the  angle  between  L\  and  Li, 
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J 


Figure  11:  Structure  of  the  Voronoi  Surface  of  Cylinder  and  Plane 


and  the  angle  between  L\  and  L 3.  It  is  not  difficult  to  see  that  the  lines  // ,  and 
7/2  intersect  in  a  point  v  on  the  axis  of  the  cylinder  g. 

Consider  any  plane  Q  through  the  cylinder  axis,  and  note  that  (7  is  not 
necessarily  perpendicular  to  /.  Then  Q  contains  v  and  a  point  w  in  which  (.’ 
intersects  Q.  Note  that  both  points  lie  on  k.  We  argue  that  the  line  177 ~  lies 
on  the  Voronoi  surface  7i.  Now  the  shortest  distance  of  a  point  g  on  this  line, 
from  the  cylinder  g,  varies  linearly  with  the  distance  of  q  from  w.  Similarly, 
the  shortest  distance  of  q  from  the  plane  /  varies  linearly  with  the  distance  of 
7  from  te,  which  is  evident  when  considering  the  plane  Q'  through  the  line  v ,  tr 
that  is  perpendicular  to  /.  See  also  Figure  11.  Since  both  u  and  w  are  on  h, 
and  since  the  distances  from  /  and  g  vary  linearly,  the  two  distance  functions 
must  be  the  same.  That  is,  the  line  v~  w  is  on  7». 

We  have  shown  that  h  contains  all  lines  through  points  of  C  and  through  t\ 
It  follows  that  h  is  a  cone  with  vertex  v  and  base  line  C. 

Case  ( i v ):  Assume  that  /  contains  the  axis  of  g,  so  g  intersects  /  in  the 
degenerate  conic  C  consisting  of  two  intersecting  lines.  Note  that  C  is  on  h. 
Consider  a  point  p  on  h  other  than  the  vertex  v  of  the  cone.  The  point  of 
closest  approach  to  p  on  /  is  pj  and  on  g  it  is  pB.  We  consider  the  line  /,/ 
through  v  and  p/,  the  line  L3  through  pg  and  v,  and  the  line  Lh  through  p 
and  v.  We  want  to  show  that  L ^  is  on  the  Voronoi  surface.  Now  Lj  is  in 
/  and  Lg  in  g.  Let  7(A)  =  Ap  +  (1  -  A)u  be  any  point  on  L an  d  7/i  A )  = 
A  pj  +  (1  -  X)v  and  7S(A)  =  A  pg  +  (l  -  A)i»  the  corresponding  points  on  I.j 
and  Lg.  Dy  proportionality,  the  distance  (7(A), 7/(A))  is  equal  to  the  distance 
(7(A),  73(A)),  and  the  angles  /(?( A), 7/(A),  v)  and  Z (g( A ) ,  ( A ),  v)  both  aie  right 

angles.  Hence,  the  line  L *  is  on  the  Voronoi  surface;  i.e.,  the  Voronoi  surface  is 
a  cone. 
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Figure  12:  Section  of  the  Voronoi  Surface  of  a  Cone  and  a  Plane 

We  now  show  that  the  curve  in  which  h  intersects  a  plane  P  perpendimini 
to  /  is  a  hyperbola.  Consider  a  point  p  on  h  at  distance  d  from  both  /  and  g. 
Then  p  lies  on  a  plane  parallel  to  /,  at  distance  d ,  that  intersects  P  in  a  line 
If.  Moreover,  p  lies  on  a  cone,  of  equal  apex  angle  a  as  g,  that  intersects  P  in  a 
circle  cg.  If  rp  is  the  radius  of  the  circle  in  which  g  intersects  P ,  then  the  radius 
of  Cg  is  rp  +  d/cos(a).  Referring  to  Figure  12,  we  have 

(rp  +  d/  cos(a))2  =  x7  +  d2 

Algebraic  manipulation  brings  this  equation  into  the  form 

(d/  cos(a)  -f  rP/  sin2(ct))2  -  i2  -  rpcot2(a)  =  0  (2) 

For  a  fixed  plane  P,  rp  and  a  are  constants,  hence  (2)  is  the  equation  of  a 
hyperbola. 

If  /  does  not  contain  the  axis  of  the  cone  g,  nor  the  vertex  of  </,  then  /  and  </ 
intersect  in  the  nondegenerate  conic  C.  We  consider  a  plane  P  perpendicular  to 
/  containing  the  axis  of  g.  P  intersects  /  in  the  line  L\,  and  g  in  the  intersecting 
lines  Li  and  L3.  Clearly,  the  angle  bisectors,  H\  of  the  lines  L\  and  Li,  and 
Hi  of  the  lines  L\  and  L3  are  on  the  Voronoi  surface  h.  The  intersection  of 
Hi  and  JI 2  is  on  the  axis  of  g,  because  their  intersection  point  v  has  equal 
minimum  distance  from  both  Li  and  L3.3  Therefore,  the  same  argument  as  in 
the  cylinder/plane  case  shows  that  h  is  a  cone  with  base  line  C. 

Finally,  if  /  contains  the  vertex  of  g  but  not  its  axis,  then  h  is  a  cone  whoso 
vertex  is  the  vertex  of  G.  To  argue  this  case,  consider  an  offset  of  /  by  some 
distance  d  and  a  plane  parallel  to  /,  at  distance  d.  Their  intersection  is  a 

3Thi»  Argument  also  Applies  to  cAse  (iii). 
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Figure  13:  Voronoi  Surface  when  Plane  /  is  perpendicular  to  Torus  axis 


nondegenerate  conic  C  that  lies  on  h  and  is  the  base  of  the  cone.  The  details 
are  straightforward. 

Case  (v);  By  Proposition  4.5,  h  is  a  quartic  surface.  Let  L  be  the  axis  of 
symmetry  of  the  torus,  and  consider  any  plane  P  containing  L.  P  intersects  <j 
in  two  circles  C\  and  C?,  and  /  in  a  line  Lf.  Let  p  be  a  point  in  P.  Since  / 
is  perpendicular  to  L,  the  point  p/  in  /  minimizing  the  distance  p,)>j  is  in 
Moreover,  since  L  is  the  axis  of  symmetry  of  g,  the  point  pg  on  g  minimizing 
the  distance  p,pg  is  also  in  P.  See  also  Figure  13.  Therefore,  P  intersects  the 
Voronoi  surface  h  in  curves  that  must  be  the  Voronoi  curves  of  C\  and  Lj,  and 
of  C 2  and  Lj.  These  curves  are  parabolas.  Since  the  distance  of  the  circles 
from  the  axis  L  is  fixed,  h  is  a  parabolic  torus  obtained  by  rotating  the  Voronoi 
curves  about  L.  □ 

Theorem  4.7 

Let  h  =  Vor (f,g),  where  /  is  a  sphere. 

(i)  If  g  is  a  sphere  of  equal  radius,  then  h  is  a  plane  or  a  prolate  spheroid.  If  g 
is  a  sphere  of  unequal  radius,  then  h  is  a  prolate  spheroid  or  a  two-sheeted 
hyperboloid. 

(ii)  If  g  is  a  cylinder  whose  axis  contains  the  center  of  the  sphere,  then  h  is 
a  plane  or  a  surface  obtained  by  rotating  a  parabola  about  a  line  perpen¬ 
dicular  to  the  axis  of  symmetry. 

(iii)  If  g  is  a  cone  whose  axis  contains  the  center  of  the  sphere,  then  h  is  a 
figure  of  revolution,  about  the  axis  of  the  cone.  The  components  of  h  arc 
obtained  by  rotating  a  parabola  around  a  line  that  is  not  perpendicular 
to  the  axis  of  symmetry.  When  g  is  tangent  to  /,  then  h  has  a  component 
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Figure  14:  Voronoi  Surface  of  Cylinder  and  Sphere 


that  is  a  right  circular  cone. 

(iv)  If  g  is  a  torus  whose  axis  of  revolution  contains  the  center  of  the  sphere, 
then  h  is  a  figure  of  revolution.  Its  components  are  obtained  by  revolving 
ellipse,  or  hyperbola  around  the  torus  axis.  If  the  radius  of  g  is  equal  to 
the  minor  radius  of  the  torus,  then  h  contains  also  a  cone,  a  cylinder,  or 
a  plane. 

Proof 

Case  (i):  Consider  a  plane  P  containing  the  centers  of  both  spheres.  If  p  is  any 
point  on  P,  then  the  points  pj  on  /  and  pg  on  g  nearest  to  p  are  also  in  /’.  It 
follows  that  h  is  a  surface  of  revolution  whose  axis  is  through  the  centers  of  / 
and  g.  The  statement  now  follows  from  Propositions  4.3  and  4.4. 

Case  (ii):  If  the  cylinder  axis  contains  the  center  of  the  sphere,  then  t lie  Voronoi 
surface  must  be  a  surface  of  revolution  about  the  cylinder  axis.  By  Proposition 

4.2,  Ii  is  obtained  by  rotating  a  parabola  about  its  latus  rectum.1  It  has  a 
spindle  as  shown  in  Figure  14.  When  the  radius  of  the  sphere  is  equal  to  the 
cylinder  radius,  by  Proposition  4.1,  the  even  Voronoi  surface  is  a  plane  through 
the  center  of  the  sphere  and  perpendicular  to  the  cylinder. 

Case  (iii):  Clearly  h  is  a  surface  of  revolution  about  the  axis  of  the  cone.  C onsid- 
cring  a  plane  P  containing  the  axis  of  g,  the  statement  follows  from  Proposition 

4.2.  See  also  Figure  15. 

Case  (iv):  Again,  h  must  be  a  surface  of  revolution  whose  geometry  is  de¬ 
termined  by  the  section  with  planes  P  containing  the  axis  of  the  torus.  By 
Propositions  4.3  and  4.4,  the  Voronoi  curves  of  the  circles  in  which  P  intersects 
/  and  g  are  ellipse,  hyperbola,  or  line.  Depending  on  the  axis  position,  the  line 
can  generate  a  cone,  a  cylinder,  or  a  plane  as  component  of  the  Voronoi  surface 
h.  □ 

*The  line  through  the  focus  of  the  parabola  perpendicular  to  the  axis  of  symmetry. 
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Figure  15:  Voronoi  Surface  of  Cone  and  Sphere 


Theorem  4.8 

Let  k  =  Vor(/,g),  where  /  is  a  cylinder. 


(i)  If  g  is  a  cylinder  of  equal  radius,  then  the  even  Voronoi  surface  is  a  hyper¬ 
bolic  paraboloid  if  the  axes  are  skew,  a  pair  of  planes  if  the  axes  intersect, 
and  a  single  plane  if  the  axes  are  parallel.  If  g  is  a  cylinder  of  different 
radius  with  a  parallel  axis,  then  h  is  an  elliptic  or  a  hyperbolic  cylinder. 

(ii)  If  g  is  a  cone  and  the  axes  of  /  and  g  coincide,  the  odd  and  even  Voronoi 
surfaces  are  cones. 

(iii)  If  g  is  a  torus  and  the  axes  of  /  and  g  coincide,  then  /i  is  a  paraboloid  of 
revolution.  If  the  cylinder  diameter  does  not  exceed  the  smaller  radius  ol 
the  torus,  then  the  Voronoi  surface  does  not  contain  a  spindle,  otherwise 
it  does. 

Proof 

Case  (i):  Assume  that  /  and  g  have  equal  radius.  If  the  axes  of  /  and  ij  are 
parallel  and  not  coincident,  then  h  is  a  plane  by  Propositions  4.3  and  4.1.  If  the 
axes  intersect,  then  the  intersection  curve  of  /  and  g  is  a  pair  of  ellipses  that  lie 
in  two  planes,  say  Pi  and  Pj.  Note  that  these  planes  bisect  the  angle  between  the 
intersecting  axes  of  /  and  g.  Consider  enlarging  the  radius  of  both  cylinders  by 
the  same  distance  d.  Evidently,  the  intersection  curve  of  the  enlarged  cylinders 
is  on  /i,  and  is  also  in  Pi  and  Pj.  Hence  h  consists  of  Pi  and  Pi. 

Assume  that  /  and  g  have  equal  radius  but  skew  axes.  Without  loss  of 
generality,  the  equation  of  /  is  x2  -f  y2  —  r2  =  0,  and  the  equation  of  <j  is 
9\{xiViz)  —  r2  =  0,  where  g\  is  a  quadratic  polynomial  in  x,y,x  without  a 
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Figure  16:  Voronoi  Surface  Point  for  two  skew  Cylinders 

constant  term.  The  even  Voronoi  surface  is  therefore  obtained  from 

x2  +  y2  -  (r  +  d)2  =  0 

gi(x,y,z)  -  (r  +  d)2  =  o 

by  elimination  of  d,  efTected  by  subtraction.  Since  x2  +  y2  -  <j\  is  quadratic, 
the  even  Voronoi  surface  is  quadratic.  Let  /  be  the  line  of  closest  approach  of 
the  axes  of  f  and  g,  and  assume  that  the  axes  of  /  and  g  pass  at  right  angle, 
at  distance  m.  The  plane  P  containing  both  /  and  the  axis  of  /  intersects  /  in 
a  pair  of  lines  and  g  in  a  circle.  It  is  not  hard  to  see  that  h  intersects  l’  in  a 
parabola.  Dy  symmetry,  h  intersects  the  plane  P'  containing  l  and  the  axis  ol 
g  also  in  a  parabola.  Next,  consider  a  plane  P\  parallel  to  P ,  at  distance  r.  Let 
p  be  a  point  on  the  intersection  of  h  and  P\,  at  distance  d  from  both  cylinders. 
Referring  to  Figure  16,  we  have  (m  4-  x)2  4-  y2  =  (r  +  d)2  and  c2  +  x2  -  (  v  4-  <l)2 ■ 

Subtracting,  we  obtain  m 2  -j-  2 mx  -f-  y2  -  c2  =  0,  hence  h  intersects  /’ |  in  a 

parabola.  It  follows  that  h  is  a  hyperbolic  paraboloid.  If  the  angle  between  the 
skew  axes  is  not  a  right  one,  then  a  purely  algebraic  argument  can  be  given, 
roughly  as  follows.  Consider  /i  =  x2  +  y2  —  g\  and  compute  its  invariants.  The 
invariants  certify  that  h  is  a  hyperbolic  paraboloid  in  this  case  also. 

If  /  and  g  have  different  radii  but  have  parallel  axes,  then  the  statement 
follows  from  Propositions  4.3  and  4.4. 

Case  (ii):  Consider  the  even  Voronoi  surface  h.  it  is  a  surface  of  revolution 
about  the  common  axis.  Consider  the  plane  P  through  the  axis.  It  intersects  / 
in  two  parallel  lin  s  L\  and  L?  and  g  in  two  intersecting  lines  I3  and  L4  that 
meet  at  a  point  v  on  the  common  axis.  The  bisectors  H\  and  II 2  of  the  angles 
formed  by  L\ L3  and  L2L4  are  iines  consisting  of  points  that  are  equidistant  from 
/  and  g.  Also  II \  and  II 2  meet  at  a  point  u  on  the  common  axis.  The  surface 


24 


of  revolution,  h,  formed  by  //t  and  II2  is  a  cone  with  vertex  at  u.  Moreover, 
since  the  two  triangles  on  P  with  vertices  on  L\  and  L 2  and  the  common  base 
11 v  are  equal  in  ail  respects,  the  apex  angles  of  h  and  g  are  also  equal.  It  is  easy 
to  see  that  the  odd  Voronoi  surface  is  also  a  cone  since  it  is  generated  bv  the 
perpendiculars  to  //t  and  //2  at  the  points  of  intersection  with  L\  and  L2. 

Case  (iii):  Consider  any  plane  P  containing  the  coincident  axes  of  /  and  g. 
The  point  p/  on  /  of  closest  approach  to  a  given  point  p  in  P  must  lie  in  P, 
and  the  same  is  true  for  the  point  pg  on  g.  Hence  h  is  a  figure  of  revolution. 
P  intersects  /  in  two  lines  parallel  to  the  axis,  and  g  in  two  circles.  If  g  and 
/  intersect,  By  Proposition  4.2,  P  intersects  h  in  a  pair  of  intersecting  and  a 
pair  of  non-intersecting  parabolas.  Since  the  lines  in  which  P  intersects  /  are 
parallel  to  the  axis  of  /,  the  parabolas  revolve  around  a  line  perpendicular  to 
the  line  of  symmetry. 

If  the  cylinder  diameter  is  less  than  the  smaller  diameter  of  the  torus,  h  is  a 
paraboloid  of  revolution.  Moreover,  the  parabola  orientation  is  such  that  there 
is  no  spindle.  If  the  cylinder  diameter  is  larger  than  the  smaller  diameter  of  the 
torus,  then  the  odd  Voronoi  surface  contains  a  spindle.  □ 

Theorem  4.9 

Let  h  =  Vor (f,g),  where  /  is  a  cone. 

(i)  If  g  is  a  cone,  and  the  axes  of  /  and  g  are  coincident,  then  the  components 
of  h  consists  of  cones. 

(ii)  If  g  is  a  torus  and  the  axes  of  /  and  g  coincide,  then  h  is  a  figure  of 
revolution  obtained  by  rotating  a  parabola  about  a  line  not  perpendicular 
to  the  axis  of  symmetry. 

Proof 

Case  ( i ):  If  /  and  g  have  the  same  axis,  then  h  is  clearly  a  figure  of  revolution. 
A  plane  P  containing  the  axis  intersects  /  and  g  in  two  pairs  of  intersecting 
lines,  and  the  bisectors  of  two  lines,  one  from  /  and  the  other  from  g,  are  the 
generators  of  h.  Therefore,  h  lias  right-circular  cones  as  components. 

Case  (ii):  Clearly  h  is  a  figure  of  revolution.  Let  P  be  a  plane  containing  the 
common  axis  of  /  and  g.  P  intersects  /  in  two  lines  and  g  in  two  circles.  By 
Proposition  4.2,  P  intersects  h  in  four  parabolas.  Since  the  lines  in  which  P 
intersects  /  are  at  an  angle  with  respect  to  the  axis  of  /,  none  of  the  parabolas 
revolves  around  a  line  perpendicular  to  the  line  of  symmetry.  □ 

Theorem  4.10  Let  h  =  Vor(/,<7),  where  /  and  g  are  tori  with  coincident  axes. 
Then  h  is  a  figure  of  revolution  obtained  by  rotating  an  ellipse  or  a  hyperbola 
about  a  line  perpendicular  to  one  of  their  axes.  If  the  minor  radii  of  the  tori 
are  equal,  then  the  even  Voronoi  surface  is  a  cylinder  or  a  cone.  If  /  and  g 
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are  congruent  with  a  common  center,  then  the  even  Voronoi  surface  is  a  pair  of 
planes  containing  the  intersection  curve. 

Proof 

Obviously  h  is  a  figure  of  rotation.  Hence  the  the  first  two  parts  of  theorem 
follows  from  Propositions  4.3  and  4.4.  If  /  and  g  are  congruent  and  have  the 
same  centers,  then  their  intersection  consists  of  two  plane  curves,  in  the  planes 
Pi  and  P-i-  The  d-oflsets  of  the  tori  are  also  congruent  with  common  centers, 
and  their  intersection  curves  lie  in  P\  and  Pj.  Hence  vhe  Voronoi  surface  lias 
two  planes  as  components.  By  symmetry,  there  is  no  additional  component.  □ 


5  Conclusion 

We  have  presented  a  conceptually  simple  method  for  constructing  the  interior 
skeleton  of  CSG  objects.  The  individual  elements  of  the  skeleton  were  shown  to 
be  the  well  known  Voronoi  curves  and  surfaces.  In  general,  the  Voronoi  elements 
have  a  complicated  geometry.  We  analyzed  a  number  of  special  cases  in  which 
the  Voronoi  surface  of  pairs  of  CSG  surfaces  has  a  simple  geometric  structure. 

Skeleton  surfaces  have  been  proposed  for  use  in  a  variety  of  important  appli¬ 
cations  such  as  computer-aided  design,  motion  planning,  mesh  generation  and 
shape  analysis.  Our  investigations,  while  illuminating  the  intrinsic  geometric 
complexity  of  such  surfaces,  provides  insights  that  can  be  exploited  in  skeleton 
construction  algorithms. 
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